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The distribution of Galactic-disk wide binaries shows a clear break in slope 
at projected separations of about r± ~ 2500 AU in two basically independent 
surveys by Chaname & Gould and Lepine Sz Bongiorno. The latter also showed 
that the frequency of wide-binary companions to G-star primaries declines mono- 
tonically as a function of companion mass. We show that both effects can be ex- 
plained by the operation of Heggie's law in the typical open-cluster environments 
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where binaries form. One immediate conclusion is that most Galactic-disk stars 



formed in open clusters with internal dispersions of a few hundred meters per 
^ . second and disruption times of a few hundred million years. 
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Two separate studies demonstrate that there is break in the distribution of Galactic-disk 
wide-binary separations at about r± ~ 2500 AU. For large separations, AO > 30", Chaname 
& Gould (2004) find a power-law distribution dN/dA0 ~ A0~ a , with a = 1.67 ±0.07, while 
for 10" < AO < 30", they find a flat "Opik's (1924) Law" distribution, i.e., a ~ 1. Because 
the characteristic distance of the Chaname & Gould (2004) sample is about d ~ 60 pc, this 
angular-scale break point corresponds to a physical scale of r± ~ dA9 ~ 1800 AU. Lepine & 
Bongiorno (2006) studied wide-binaries with Hipparcos (ESA 1997) primaries. Because their 
targets have parallaxes, they present their results directly in terms of r± (rather than A6). 
They find a break at about r± ~ 3000 AU. Each estimate is uncertain at about the 20% 
level. Lepine & Bongiorno (2006) are hampered by small number statistics near the break 
because their bright-primary sample does not permit them to probe inward of A9 < 20". On 
the other hand, while Chaname & Gould (2004) probe to smaller angular separations and 
have much better statistics, it is difficult to make precise the interpretation of their angular 
break in terms of physical scales because of the substantial uncertainties in the individual 
distances of their sample. Thus, the two studies are broadly consistent. Moreover, as we will 
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show below, one actually expects the Chaname & Gould (2004) sample to have a somewhat 
smaller break point. 

Neither of these two surveys can probe inward of about r± = 1000 AU. However, it 
is known from the work of Duquennoy & Mayor (1991) that the G-star binary distribution 
roughly obeys Opik's Law over about 4.5 decades of separation, from 10~ L5 AU to 10 3 AU. 
(The Duquennoy & Mayor 1991 distribution is often expressed as a Gaussian, but it should 
be remembered that the abscissa of this "Gaussian" is in log while the ordinate is linear.) 

For binaries observed at random orientations, their semi-major axis a is related to 
their projected separation by (cr 1 ) = (7r/4)(r^ 1 ), so the observed break at r± ~ 2500AU 
corresponds to a ~ 3000 AU, i.e., to orbital velocities t> or b ~ 0.5 km s -1 (assuming a typical 
total mass M ~ 1 M ). 

Lepine & Bongiorno (2006) suggested that the observed break could be due to dynamical 
effects in the clusters in which the binaries were born. Indeed, the orbital velocity at abreak ~ 
3000 AU is very similar to the 1-dimensional dispersion observed in open clusters, which 
typically ranges from o ~ 0.3 km s -1 to a ~ 1 kms -1 . Hence, it immediately argues that the 
break is an effect of "Heggie's (1975) Law", which states that "hard binaries get harder, soft 
binaries get softer" . The boundary point of this law is defined by the internal binding energy 
of the binary being equal to the mean kinetic energy of the ambient perturbers. Since, the 
binary components and the perturbers have roughly the same mass, the break point occurs 
when t> orb ~ a. 

Consider a "soft binary" population whose members are being injected with energy at 
a rate that is a function of their binding energy E^. dE^/dt ~ E^ '. If the population 
is subjected to this process for a sufficiently long time, it will reach a steady state with 
dN/dE b ~ £f, corresponding to a separation distribution dN/da ~ a" 2_/3 . That is, a = 
2 + (3. To a good approximation, the injection of energy into the binary is independent 
of separation, i.e., f3 = 0. This would predict a = 2, somewhat larger than the value 
a = 1.67 ± 0.07 observed by Chaname & Gould (2004). To next order, the energy injection 
actually grows logarithmically with separation, so (3 > 0, which goes in the wrong direction 
but only by a small amount. Hence, if the argument we are giving is correct, the explanation 
for the discrepancy in slopes must be that the cluster dissolves before it has time to reach 
its "asymptotic state" . 

The timescale for binary evolution is T ~ (nSf) -1 , where n = p/m is the ambient 
perturber number density, p is the mass density, m is the mass of a typical perturber, 
£ = ^brcak is the cross section for a major perturbation at the break point ab rC ak ~ 3000 AU, 
and v — a is the ambient velocity. Using the virial theorem, AirGp ~ (a/R) 2 , one finds 
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T ~ R 2 j (abreak^) ~ 100 Myr, where we have assumed M ~ 2m, adopted a cluster radius of 
R ~ 1 pc, imposed Heggie's Law (i> orb ~ a), and dropped factors of order unity. Since open 
clusters generally dissolve on timescales that are one or several times this binary-disruption 
timescale, it is plausible that the binary distribution does not have time to fully reach its 
asymptotic state. 

The explanation just given makes an important prediction. The binding energy of a 
binary scales Ej, ~ MiM 2 /a, where Mi and M 2 are the component masses. The Hipparcos 
primaries in the Lepine & Bongiorno (2006) study are virtually all solar- type stars, i.e., 
Mi ~ Mq. This means that for secondaries of different masses, the break point scales as 
Obrcak ~ M 2 . It is plausible to assume that the secondaries with a ~ abreak are initially 
drawn randomly from the field population. Then we would predict that after the binaries 
diffuse to larger a, the ratio of secondaries to field stars of the same mass would fall by 
a factor a^ Teak , or in other words as M 2 . Figure 9 of Lepine & Bongiorno (2006) shows 
the frequency of observed secondaries compared to what would be expected based on a 
distribution normalized at 4 < My < 8, i.e., stars of mass M 2 = 0.8 Mq. The foregoing 
argument would predict that at My = 12 (M 2 ~ 0.25 M©), the observed secondaries should 
be deficient by a factor ~ (0.25/0.8) 1 ' 67 ~ 0.15. The actual deficiency is about 0.5, which 
is significantly less dramatic. Nevertheless, this figure does show the expected overall trend. 
One possible explanation for the discrepancy is that the timescale for disruption of the 
binaries with smaller secondaries at their break point is considerably longer simply because 
their orbits present smaller cross sections. Since the binary disruption timescales at "typical" 
masses are already of order the cluster-disruption timescale, this increase in binary-disruption 
timescale could substantially mitigate the accelerated disruption relative to the naive scaling 
we have given. 

The same argument predicts that the Chaname & Gould (2004) sample should have 
a smaller ab rea k than the Lepine & Bongiorno (2006) sample because each of the latter is 
guaranteed to have a Hipparcos (i.e., roughly solar mass) component, and so to have a 
systematically higher binding energy. 

We attempt a first test of our hypothesis by dividing the sample shown in Lepine & 
Bongiorno (2006) into three subsamples, with M v < 8, 8 < M v < 12, and M v > 12. See 
Figure 1. Our prediction would be that the fainter stars should plateau at smaller r±. If 
there is any trend it would appear to go in the opposite direction, although these subdivided 
data are quite noisy and could still be subject to selection effects if the faintest stars are 
more difficult to detect at ~ 20" separations that Lepine & Bongiorno (2006) believe. 

An important implication of this argument is that the break in the binary separation 
function tells us about the typical conditions in which disk stars form, namely that over 
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the lifetime of the Galactic disk, most stars formed in clusters with velocity dispersions of 
order a few hundred meters per second and that these clusters disrupt on timescales of order 
a few hundred million years. Figure 10 of Chaname & Gould (2004) hints that the halo 
binary breakpoint occurs at at least a somewhat smaller semi-major axis than that of the 
disk. If this is confirmed by future observations, it would imply that halo stars were born 
in environments that were kinematically at least slightly hotter than disk stars. (Note that 
it is also possible, in principle, to produce a smaller ab re ak by selecting a sample of halo 
binaries with systematically smaller companion masses. However, the Chaname & Gould 
2004 sample actually has the opposite bias: since the halo stars are a factor ~ 4 farther 
away, they are actually biased toward being more luminous - and so more massive - stars.) 

Detailed verification of the scenario presented here will require simulations that simul- 
taneously model the disruption of the binary and of the cluster in which it initially resides. If 
verified, the same simulations will permit more precise characterization of the typical clusters 
in which today's field binaries were born. 

We thank Scott Gaudi for valuable comments on the manuscript. This work was sup- 
ported by NSF grant AST 042758. Any opinions, findings, and conclusions or recommenda- 
tions expressed in this material are those of the authors and do not necessarily reflect the 
views of the NSF. 
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Fig. 1. — Frequency of binary companions as a function of projected separation in three 
luminosity bins for the secondaries: My < 8 (triangles), 8 < My < 12 (squares), and 
My > 12 (stars). These correspond roughly to secondary mass intervals M2 > 0.65 M , 
0.25 M < M 2 < 0.65 M , and M 2 < 0.25 M . For clarity, the 3 sets of points are slightly 
offset in the horizontal direction. The data are taken form Lepine & Bongiorno (2006) and 
correspond to their Figure 10. According to the argument presented here, the lower-mass 
bins should plateau at smaller projected separations. The data do not show this trend but 
they are noisy and may still be affected by selection. 



